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Abstract 



We present an approximative calculation of the ground-state en- 
ergy for the anisotropic oscillator with a potential 



Using an instanton solution of the isotropic action 5 = 0, we obtain 
the imaginary part of the ground-state energy for small negative g as 
a series expansion in the anisotropy parameter 5. From this, the large- 
order behavior of the ^-expansions accompanying each power of 5 are 
obtained by means of a dispersion relation in g. These g-expansions 
are summed by a Borel transformation, yielding an approximation to 
the ground-state energy for the region near the isotropic limit. This 
approximation is found to be excellent in a rather wide region of 6 
around 5 = 0. 

Special attention is devoted to the immediate vicinity of the 
isotropic point. Using a simple model integral we show that the large- 
order behavior of an 5-dependent series expansion in g undergoes a 
crossover from an isotropic to an anisotropic regime as the order k of 
the expansion coefficients passes the value fe cr0 ss ~ VH- 



V(x, y) = \{x 2 + y 2 ) + I [x 4 + 2(1 - 5)x 2 y 2 + y 4 
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1 Introduction 



Phase transitions in anisotropic systems with cubic symmetry have attracted 
much interest in the literature [0-0. Especially well studied are correspond- 
ing models in quantum mechanics. To gain an analytic insight into the latter, 
Banks, Bender, and Wu (BBW) studied a Hamiltonian with a potential 



Using multidimensional WKB techniques they derived the large-order behav- 
ior of the perturbation series for the ground-state energy 



function of the anisotropy parameter 5. 
In 1990, Janke derived the same results with more efficiency from a 
path integral for the imaginary part of the energy E. The imaginary part 
contains information on the tunneling decay rate of the ground state for 
g < 0, and determines directly the large-order behavior of the perturbation 
coefficients via a dispersion relation in the complex coupling constant plane. 
Both, BBW and Janke, find a different large-order behavior of the isotropic 
system 5 = and the anisotropic system 5^0. They do not discuss, 
however, the interesting question of how the latter goes over into the former 
as 5 goes to zero. 

It is the purpose of this paper to fill this gap. For an optimal understand- 
ing of the expected behavior we shall not attack directly the path integral 
involving the potential ([!]), but first only the corresponding simple integral. 
For this integral, a perturbation expansion of the form (|2]) yields exactly- 
determined 5-dependent perturbation coefficients. The coeffients E k (5) are 
shown to have a large-order behavior which undergoes a crossover between 
the earlier derived isotropic and anisotropic behaviors when the order k passes 
the crossover value k CVQSS ~ 

The expansion terms of a model integral with the potential ([!]) counts 
the number of terms in a perturbation expansion of the quantum mechanic 
and the field theory. Thus the bare model integral is sufficient to derive 
nontrivial information on the large-order behavior of the eventual object of 
interest, quantum field theory. It turns out that for resumming the (/-series, 



V= | [x 4 + 2(l-5)xV + / . 



(1) 



E = J2E k (5)g k 



(2) 
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asymptotic large-order estimates for the 8- dependent coefficients can be used 
only in the anisotropic regime k\8\ 3> 1. In the isotropic regime k\8\ <C 1, 
on the other hand, it is impossible to truncate the large-order expansion 
of the perturbation coefficients after a finite number of terms. Thus the 
neighbourhood of the isotropic system 5 = needs an extra investigation. In 
the context of quantum field theory, this was recently delivered in ||. 

The imaginary parts of physical quantities at small negative g can be 
calculated with the help of classical solutions called instantons. In systems 
sufficiently close to the isotropic point it is not necessary to know the exact 
instanton solutions for all 5. The knowledge of the solution at the symmetry 
point 5 = is perfectly sufficient, around which the imaginary parts can be 
expanded in powers of 8. 

After having understood the model integral, we shall perform the same 
analysis for an anisotropic quantum mechanical system, which represents an 
one-dimensional 4 -field theory with cubic anisotropy. 

The paper is organized as follows. In Section 2 we develop a simple 
resummation procedure by which the divergent power series expansion of 
a function Z(g) = J2k^kg k is converted into an almost convergent series 
Y^ p a>plp{g)- Here I p {g) are certain confluent hypergeometric functions which 
possess power series expansions in g with similar large-order behavior as 
the system under study. In Section 3 we shall analyse the above-mentioned 
crossover in the large-order behavior for the simple model integral. In par- 
ticular, we shall justify the resummation procedure of Section 2 and the 
methods in H] to be the perfect tools in approximating the integral for the 
region near the isotropic limit 8 — > 0. In Section 4, finally, we present a 
similar calculation for the ground state energy of the anharmonic potential 
with cubic anisotropy. 

In addition to this more standard resummation procedure we analyse the 
model also within the variational perturbation theory developed in 



Variational perturbation theory yields uniformily and exponentially fast con- 
verging expansion for quantum mechanical systems with quartic potentials 
[|T3]]. The uniform convergence was first proven for the partition function 
of the anharmonic integral, later for the quantum mechanical anharmonic 



oscillator with coupling strength g in several papers [ 14J . Recently, the proof 



was sharpened and extended to the energies |T3 



The input for the quantum mechanical model is provided by the ex- 
act Rayleigh-Schrodinger perturbation coefficients of the ground-state en- 
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ergy which we derive from an extension of recursion relations first shown by 
Bender and Wu (BW) pf. 

2 Resummation 

We begin with developing a practical algorithm for a Borel resummation of 
a divergent perturbation series 

Z(g) = J2Z k g k . (3) 

k 

Our method will be most efficient under the following conditions: 

1. From low-order perturbation theory we know the expansion coefficients 
Zf~ up to a certain finite order N. 

2. From semi classical methods we are in the possession of the high-order 
information in the form 

Z k k -=^ 7 (_l) fe A;!fcV(l + ^ + | + ---). (4) 

3. By some scaling arguments we are able to assure a power behavior in 
the strong coupling limit 

Z(g) ^ ng a (5) 
The idea of the algorithm is the following: 

It must be possible to construct an infinite, complete set of Borel summable 
functions I p {g) which satisfy the high-order and strong-coupling conditions 
(J3j) and (|5]). These functions can be used as a new basis in which to reexpand 

Z{g): 

oo 

z (g) = J2 a P I p(9) , (6) 

The series (^) should be such that the knowledge of the first (N + 1) coeffi- 
cients in the power series expansion (|]) is sufficient to determine directly the 
first (N + 1) coefficients a p , yielding an approximation 

N 

Z( g )*Z( N \g) = j:a p I p (g)- (7) 

p=0 
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This would then be a new representation of the function Z(g) with the same 
power series up to g N but which makes use of large-order and strong-coupling 
informations (||) and (0). In the limit of large N, the series (|7|) is expected 
to converge towards the exact solution. 

The functions I p {g) being Bore summable have a Borel representation 

I p (g)= / dte-H b0 Bl°(gt). (8) 

J 

Parametrized by some b$ and integer p, what are the conditions on B p (gt), 
such that I p (g) satisfies (|J) and fl5|) for all p? The answer is most easily found 
with the help of the hypergeometric functions 

2 F 1 (a, b; c; -agt) = ^ — r, (9) 

k=0 

with appropriate parameters a(p), b(p) apd c(p). The Pochhammer symbol 
(a)k stands short for (a)k = T(a+k) /T(a). These functions have the following 
virtues: First, they are standard special functions of mathematical physics 
whose properties are well-known. Second, they have a cut running from 
t = —l/\ag\ to minus infinity which is necessary to generate the large-order 
behavior Third, they have enough free parameters to fit all input-data. 
The first property permits an immediate calculation of the Borel integral @, 
which is simply a Laplace transformation of t b ° 2 i<i (a, b; c; —cgt) 

roc T(c) 

/ dte-H b0 2 F 1 {a, b- c- -agt) = —U—E(a, b,b + l:c: l/ag) (10) 
Jo 1 (a)L (b) 

The resulting E(a,b,b + 1 : c : I /erg) is MacRobert's .E-function. Using its 



asymptotic expansion (see Ref. [|T7[], page 203) it is easy to verify that our 
ansatz reproduces the large-order behavior (|). Indeed, for large k the power 
series 

T(r) 00 
—AJ-E(a, b,b + l:c: l/ag) ee £ e k g k (11) 

r(a)r(6) ^ 

has coefficients which grow like 

e k ^ yW^) k] k a+b - c+b ^a k . (12) 
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Moreover, this property is unchanged if the original hypergeometric function 
is multiplied by a power (agt) p . A possible set of Borel functions are therefore 
the following functions: 

B p (gt) = (agty 2 F l {a,b;c; -agt) . (13) 



Looking at fll2] ) we see that the functions (|IBD are not completely fixed by 
a given large-order behavior. The parameter (3 in (|j) merely imposes the 
following relation upon the parameters a,b,c and b 

a + b-c + b -l=f3, (14) 

and there are many different ways to satisfy this. The specific choice will be 
suggested by practical considerations. One such consideration is that the J p 's 
should possess a simple integral representation in order to avoid complicated 
numerical work. In addition, we would like to work with parameters a, b 
and c, for which the hypergeometric function 2 F\ reduces to simple algebraic 
functions. This happens only for special sets of the parameters. A simple 
possibility is for instance (see Ref. |1§, page 556) 

1 / / \ -2a 

2 F 1 (a,a+-;2a+l;-z)=4 a [l + VTT^) (15) 

which arises by choosing the parameters a, b and c which are related by 

a + b-c = -^, ; c-26 = 0. (16) 

With this, the relation flT4] ) can be satisfied for an arbitrary value of the 
parameter a by choosing 

bo = P+\, (17) 

and we are left with only one parameter degree of freedom. This freedom 
may be used to accommodate the strong-coupling behavior of Z(g) if it is 
known. The equation ([?]) yields the condition I p (g) — > const, x g a on the 
functions I p (g)- From (H) we see that such a power behavior emerges if all 
Borel functions B p satisfy B p (z) — > const, x z a and thus 2Fi(a,b;c; —z) — > 
const, x z~ p+a [see (([3D]. The explicit representation fllSl) shows that the 
parameter a has to be taken as 

a = p — a . (18) 
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Thus we obtain the approximation Z^ N ' = Y,p=o a p Ip with 

r oo p-tfbo (rr n i\P / 







oo e~H b ° (crgt) p ( 1 
dtr^-r, — — - — 2-^1 [P ~ a iP ~ a + 7^; 2(p — a) + 1; — <7#t 



r(6o + l) 4p 



OO g 

dt 



fho 



r(6 + l)V2 2 



+ -Jl + ar?t 



2a 



2" 



1 + y/T+agi 



2p 



(19) 



where the Borel parameter bo is fixed by flTTp. The normalization constant 
l/4 p r(6o + 1) in front of the expansion functions was introduced for conve- 
nience. 

Let us now derive equations for the expansion coefficients a p in terms of 
the perturbation coefficients Z\.. All one has to do is take the asymptotic 
expansions 

OO 

ys) = E4V, (20) 

k=0 

insert these into (|7|), collect terms of equal power g k , and compare these with 
the perturbation series @. This gives the (N + 1) algebraic equations 



A? 



p=0 



k = 0,1, 



N 



(21) 



By assumption, the series on the left hand side contains only the coefficients 
a p with p < N. Thus the a p s can be computed, in principle, by inverting 
the (N + 1) x (N + 1) matrix (I)k P = Even though this can be done 
recursively for any given case, it is preferable to find an explicit algebraic 
solution for a p in terms of Z^. This is possible using the following trick. We 
rewrite the asymptotic expansion of Z^ N ' in Borel form 



N oo poo 

zW(g) = j:a p I p (g) = E4 N) 9 k = / dte~^ £ 

p=0 k=0 J0 k=0 



AN) 



(gtf 



% T(k + b + l) 



(22) 



insert the expression (]TQ|) for I p (g), and compare directly both integrands 



N 

E« 

p=0 



(I + |yi + <7gt 
' r(6 + 1) 



2a 



agt 



1 + ^/TT^gi 



^r(* + 6d + i) 



• (23) 
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Introducing the new variable 



u,ss- ^ 2 = (24) 

U _l / i i ^ 2 Jl + apt + 1 V 7 



(l + VI + <70t) VI + ^ 
we obtain from (1231) the relation valid for all a 



5^£(^W <n^- (25) 

In order to compare equal powers in w we expand on the right hand side 

(1 _ w )-W-«) = g ^ " 2 (* ~ «) ^ (26) 

which gives after a shift of the summation index from Z to p = k + I the first 
(AT + 1) coefficients a p in terms of the perturbation coefficients 



a 



p 



k=0 



(&o + l)* W \ P-fc 



(recall that = Zk for = 0, 1, . . . , N). Finally, rewriting the binomial 
coefficients by means of the identity 

{ X p) =i - 1)P { P ~ X p~ 1 )> (28) 
we obtain the more convenient expression 

p Z k (A\ k ( p + k-l-2a x 



a 



v 



k=0 



r ,: ■ pa) 



(6q + 1)* \crj \ p-k 



Thus, we have solved the original matrix inversion problem (|2T|) by translat- 
ing it to a simple problem in function theory, namely that of inverting the 
function w(agt) in fl24"|). For the purpose of calculating the integrals I p (g) 
numerically, we may use the variable w itself as a variable of integration, and 
rewrite the integral representation for I p (g) in the form: 

(±\ bo+1 f 1 , (l + w)w b «+'P [ Aw 

p{9} ~\~g) Jo r(6 + i)(i- w ) 2b o+2-+3 ex P [ (i- w )2 ag 

(30) 

Together with the explicit formula (|29|) for the coefficients a p we thus have 
solved the resummation, and it is now straightforward to calculate the ap- 
proximation (0). 
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3 Model integral 



In order to set up an approximation method for an anisotropic model in the 
neighborhood of the isotropic point 5 = 0, it is instructive to study first a 
simple toy model whose partition function is defined by a two-dimensional 
integral: 



1 

2^ 



dxdyexp j — [x 2 + y 2 ^j — j x A + 2(1 — 5)x 2 y 2 + y 



(31) 

This can be interpreted as a partition function of a 4 -theory in zero 
spacetime dimensions with cubic anisotropy. Introducing polar coordinates 
x = rcos<£> and y = rsimp, we obtain the more convenient form of the 
integral ([£!]): 



Z = — , 

2tt Jo 



OO /■27T 



dpdipexp —p — G(g,5,(p)p" 



(32) 



with p = r 2 /2 and 



G(g,S,if) = g 



6 -sm 2 (2tp) 



(33) 



After an integration over the angle (p, we find the integral 

. (5 . 



dpexp 



6 
4 



P 



4 



gp 



(34) 



where Iq(x) is a modified Besselfunction I v (x) for z/ = 0. Eq.(p4]) is useful 
for a numerical calculation of Z(g,S). It will serve as a testing ground for 
our approximations. 

Thanks to the special spacetime dimensionality of the model, the pertur- 
bation expansion of Z(g, 5) can be obtained explicitely and we can calcu- 
late the large-order behavior without doing the saddle point approximation, 
which is unavoidable in quantum mechanics and field theory. 

At first glance it seems useful to expand Z(g,S) in the form 



Z(g,6) = J2Z k (S)g h 

k=0 



(35) 
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where the perturbation coefficients are parametrized by the anisotropy 8. 
The coefficients Z k (8) may be found by expanding the integrand of fl32| ) in 
powers of g and performing the integral term by term: 



Z k (5) 



k\ 



-r(2fc + i) 1 



k/2 



Pi, 



' 4-8 N 
2 v / 4^25, 



(36) 



where P k (x) are the Legendrepolynomials. 

In Figure p] we have plotted the order dependence of these coefficients for 
the anisotropy parameter 8 = 10~ 2 . What we can see is a crossover of the 
large-order behavior from an isotropic to an anisotropic regime in the vicinity 
of an special crossover value fc croS s ~ Vl^l = m ^ ne anisotropic regime 
k\8\ ^> 1, the large-order parameter j3 has the value j3 — — 1. For k\8\ <C 1, 
however, we can read off the large-order behavior of the isotropic case, i. e. 
(3 = —1/2 (see also Figure 0). 

Using the duplication formula for Gamma-functions and the expansion 

T(k + e + l) = k\k £ [l + 0{l/k)}, (37) 
the large-/c behavior of T(2k + 1) is given by 

T(2A; + 1) = ^ l ' 2 4 k (k\) 2 k- 1 ' 2 [l + 0{l/k)\ . (38) 
An approximation of the Legendre Polynomials Pk(x) for large k including 



contributions of the order 0(1/ k) can be derived from Hobson (see Ref. JLS 
page 305): 



P k (x) = (2itk)- 1 ' 2 (x 2 - l) 1/4 (x + Vx^l) 



k+l/2 



X 



1 



\[x 1 ~ 



1 (x + \fx 2 



Substituting 

we obtain for 8 > 
„ / 4-8 N 



8ky/x 2 - 1 (x + V 
4-8 



x- 



+ 0(l/k 2 



x 



2V4^28 



(39) 



(40) 



2y/4 
1 r 2 

71 



28 J 

/ir l/2 r l/2 



k/2 



l-8/2 / 



1 + 



1 f4 -38" 



k8 



O 



( 1 



\k 2 8 2 . 



(41) 
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Figure 1: 
in Eq. 

09 = - 
5 = 10 
(In /c cros 



Crossover of large-order behavior of the expansion cofficients 
from the isotropic regime {(3 = —1/2) to the anisotropic regime 

Zk/(— 4) k k\ for the anisotropy 
~ l/\6\ = 10 2 



1). Plotted is the function f(k) = In 
~ 2 . In this case the crossover value is given by k, 
,»4.6). 
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The combination of ( |38| ) and ([IT]) yields the large-order behavior of the per- 
turbation coefficients Z k (S): 



Z k {5) 



2 l/2 



7T 



1 

kS 



1 



+ 0(S) 



+ o 



( 1 



(42) 



A similar calculation can be done for 5 < with the result: 

(2 - 5) 1 ' 2 



Z k {5) 



71 



i-l) k (A - 25) k k\k-\-5)- 1/2 



1 



+ C((5) 



+ o 



k 2 fr 



(43) 



For resumming the series fl35"D , the perturbation coefficients ( fi"2"D and (fl3|) 
can be used only for fc|<5| ^> 1. In the regime k\b \ <C 1, on the other hand, it 
is unpossible to truncate the series in ((42"|) and in fl43Tj after a finite order of 
^> 1. Thus, the isotropic regime cannot be described by resumming 
the perturbation series (|3~5p using the asymptotic results ( |42"D and (|43|). Be- 
ing interested in the region close to the isotropic limit, we therefore use an 
expansion different from fl35|), and rewrite Z(g,5) as 



oo k 

Z(g,5) = Y,Y, z *ng k 5 n 

k=0n=0 



(44) 



Then, for the reason given in ||, reasonable results should be obtained by 
resumming the g-series accompanying each power 5 n . The explicit form of 
the coefficients Z kn is given by 



J kn 



7T 



-1/2 



("I) 



1)k+n r(n + i)r(2fc + i) 

2 n r(n + l) 2 r(A;-n + l) 



(45) 



and Zkn = for k < n. For these coefficients, the expansion of the Gamma- 
functions yields the behavior for large k ^> n: 



J kn 



r (n + \) 
' 2«r(n + l) 2 



-r 



,4* 



7T 



k\k n ~ 1 ' 2 [l + 0(l/k)] 



(46) 



In the following we shall calculate the coefficients ([fj|) by means of the 
steepest descent method using the saddle points of the limit 5 — > 0. This 
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will serve as a simple preparation for the analogous method in quantum 
mechanics and field theory. 

As a function of a 6 and a complex coupling g, the integral (^) is defined 
in the half-plane Re g > 0. For Re g < 0, the integral can be calculated by 
an analytical continuation from the right into the left half-plane, keeping the 
integrand in (|31|) real. This analytical continuation can be achieved by a 
joint rotation in the complex (/-plane and of the integration contour in the 
r = (x, y)-plane. The convergence of the integral is maintained by the sub- 
stitution g — > gexp(i$) and f — > rexp(— i6/A), where 9 is the rotation angle 
in the complex g-plane. Let us assume that the function Z(g, S) is analytic 
in the g-plane, with a cut along the negative g-axis, and a discontinuity for 
g < 0. Then the rotation in the complex g-plane by an angle 9 = ±7r yields 
on the lower lip of the cut: 

Z±{-\g\,&) = z(\g\e hi ",6) . (47) 



The corresponding rotated integration contours (r^) are drawn in Figure ||[ 
The discontinuity across the cut is given by 



/■27T (j^tn f- 

Disc Z = / — / rdr exp 
jo 2n Jr 



1 4 

r n( x \ r 
-2"-%^)j 



(48) 



where the combined contour r = T + — T_ runs for r > entirely through 
the right half plane. 

In a perturbatively expansion in powers of 6, the discontinuity can be 
computed from an expansion around the saddle point 



ro = \ T1 (49) 
V \9\ 

of the isotropic case 5 = 0. Since r > 0, only the positive square root 
contributes, the negative one is automatically taken into account by the 
integration over the angle (p. Now, the contour of integration T in the right 
half-plane can be deformed to run vertically across the saddle point (see 
Figure |j), i. e. , we can integrate along a straight line: 

- • (50) 
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® 






— 7T ^ 





->- 



71 




a) 



6) 



Figure 3: Analytic continuation g — > \g\ exp(±7r): a) Rotation by angles 
±7r in the cut complex g-plane. b) Two rotated paths of integration in the 
r-plane (r > 0). 



V 



V 



\/ 



r 



\/ 



Figure 4: Deformation of the contours of integration to make them pass 
through the saddle point r . 
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The exponent in (|4*8| ) plays the role of an action, and the deviations £ may 
be considered as fluctuations around the extremal solution. The angle if is 
analogous to a collective coordinate along the motion of the instanton in the 
isotropic limit. Expanding the action up to the second order in £ around the 
extremum of the isotropic action, we obtain 



DiscZ = — - ( r-r 



x exp 



1/2 



oo p2n 



-oo JO 



m sis 



d^dip 



sin 2 (2<p) - £ 2 + O 



(51) 



Integrating out the fluctuations £ and the collective coordinate ip, and using 
the equation 

(52) 



Im Z = — Disc Z 
2i 

we obtain the following imaginary part for Z: 

/ 1 \ n+1/2 



Im Z 



n=0 



2T(n+ l) 2 \4|^| 



exp 



[i + ofo)]- 



(53) 

Each power 5 n has its own n-dependent imaginary part. Given such an 
expansion, the large-order estimates for the coefficients Z^ n (with k *C n) 
follows from a dispersion relation in g (see for example Eq. (6) in Ref. 0) 



J kn 



7T J-oo 



o , lmZ (n Hg + iQ) 
dg- 



9 



k+l 



(54) 



where Z^(g) is the coefficient of 5 n . In general if a real analytic function 
F(g) has on top of the cut along g G (—oo, 0) an imaginary part 



(3+1 



ImF(# + iO) 



-7T7 



exp 



[1 + OO/)] 



(55) 



then a dispersion relation of the form fl5"4]) leads to the asymptotic behavior 

F fc = T (-l) VA^A;! [1 + 0(1/Jfe)] . (56) 



With <j = 4 and (3 = n — 1/2, we obtain again the result fl4"6). 
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Thus, the steepest descent method using the isotropic saddle point is a 
perfect tool for calculating the large-orde behavior of the expansion coeffi- 
cients Z kn in the expansion (f£|). A great advantage of this method with 
respect to the exact calculation ( [45]) is the fact that it can be generalized 
to quantum mechanics and field theory where exact calculations would be 
unpossible. 

Before applying the resummation algorithm of the previous section we 
have to study the strong-coupling behavior, i. e., the limit of large g. This 
can simply be done by rescaling the integral ( [32"1) 



\fG(g,5,<p) 



- y 



(57) 



with G from (^) and y = pyG. Taking the limit of large g (i. e. large G) 
and integrating out the angle ip we find 



with 



Z(g,5)) 9 ^ K {S)g-W 



vr 1 / 2 00 (2nV 2 
1 ; 2 to N) 4 2 5 " 



(58) 
(59) 



Now, a resummation of the (/-series in flli]) yields a generalization of ([?]): 

N / N 



N / N \ 

Z^ N \g,5) = Y: (J2a pn l pn (g))5 n 

n=0 \p=n J 



(60) 



with the complete set of Borel summable functions 
Ipn{g) = 



eg) 



b (n)+l 



dw 



[1 + w )w bo{n)+p 



T [6o(n) + 1] (1 - W )2bo(n)+2a+3 



exp 



Aw 



w) 2 ag 



and the coefficients 



J kn 



"pn 



£ (bo(n) + l) k 



a 



p + k 



2a 



p — k 



(61) 



(62) 
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where the perturbation coefficients Z^ n are given by The parameters 

bo(n), a and a follow from the large-order behavior (|46| ) and the strong- 
coupling expansion (|5%D, respectively: 

b (n) = n + 1 
a = 4 

a = -l. (63) 

From ( |62| ) it is possible to derive the following closed formula for the coeffi- 
cients a pn : 

^- {n + 2) (-D-^tf^ ( 2 ; + : ) . (64) 

7r 1 [2n + 2) p\n\ \ p-n ) 

Inserting the exact strong-coupling parameter a = —1/2, we obtain 

r 1 (ra + l) (2ra)! 
a pn = 8»(2n+l) (n!) 2 tor ^- n ' (65) 
[ else . 

Thus we have a general result that the approximants J2p= n a pn Ipn{g) in (|60|) 
posses no terms with p > n, where n is the power of 5. In such a way the 
n-dependent functions of g associated to each 5 n are recovered exactly 

The approximation Z^ N \g,S) may then be compared with the numeri- 
cally calculated integral (p4] ) . In Figures |5] and |6| we have shown the result 
for various N and coupling constants g/4. 



4 Quantum mechanics 

After the integral model, the simplest nontrivial example for our approx- 
imation method is a 4 -theory in one spacetime dimension with a cubic 
anisotropy, which is equivalent to the quantum mechanics of an anisotropic 
anharmonic oscillator. 
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-3-2-10 1 2 

6 

Figure 5: Partition function Z of the simple integral model as a func- 
tion of the anisotropy parameter 8 with the coupling constant g/4 — 0.25. 
Comparison is made between the precise numerical result and the resummed 
perturbation series [see Eq. (|5DD] for various order N. 
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4.1 Recursions relations for the ground-state pertur- 
bation coefficients 



Consider an anharmonic oscillator with cubic anisotropy and a Hamilton 
operator 



H = - 



1 / d 2 



d 2 



2 \dx 2 dy 2 
Introducing reduced variables by a rescaling 

T 



+ Y + V 2 ) + f K + 2(1 - S)x 2 y 2 + y A ] . (66) 



x 



9 



cu 3 g 



V - 
E (o) 



-y 



(67) 



yields the dimensionless time-independent Schrodinger equation 



2 



+ 



d 2 



2 \dx 2 dy 



1 2 . 9 ( 
—r H — r 
9 A \ 



2Sx 2 y 2 



(68) 



with the associated boundary condition 

\^ n (x,y)\ — >0 , r = \jx 2 +y 2 



oo . 



(69) 



The boundery condition selects only the discrete energy eigenvalues E^ n \ 
We now consider the ground-state energy E^ = E, whose perturbation 
expansion has the form 



E 



oo oo 

EE 

m=0 l=rn 



{25) m E t 



in i 



(70) 



with the unperturbed ground-state energy E Q0 = 1. In the following, we 
refer to (|70|) as the Ray leigh- Schrodinger series and the E im as a Rayleigh- 
Schrodinger coefficient. 

In general, the ground-state energy is available from the sum of all con- 
nected Feynman diagrams having no externel legs. For an efficient compu- 
tation of the Rayleigh- Schrodinger coefficients at large orders we shall derive 
recursions for the Ei m following a method introduced by Bender and Wu 
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IHj. In this way we obtain a difference equation generating the Rayleigh- 



Schrodinger coefficients. 

Separating out the unperturbed ground-state wave function, *f? (x,y) 
exp[— (x 2 + y 2 )/2], we substitute 



*M = EE 1 (26) n exp[-(x 2 + y 2 )/2}^ kn (x,y), (71) 

n=0k=n ^ 4/ 

where &kn{x, y) is polynomial in x, y with $00 = 1- Inserting the perturbation 
expansions ([TO]) , ( ffl| ) into the differential equation (|68D, and collecting powers 
of g and 5, we find 

1 / d 2 d 2 \ ( d d\ 

2 y aor a?r J \ ox ay j 

k n k 

J2(-l) l ElO®k-l,n+J2 T,(- 1 ) lE ^k-l,n- m (72) 
1=1 m=l l=m 

Finally, the ansatz 

2fc-n 

$kn(x,y)= £ ■\;--> :2 '<f J (73) 

i,j=0 

with 

A 1 *™ = for i,j > 2k -n ; i,j,k,n <0 ; < n (74) 
gives the desired difference equation 

I kn _ (n- 1 1 \ / • 1 i \ d fcn i (n ■ , iU ■ , i \ j fcn 



2(< + j)A$ = (2i + l)(i + + (2j + l)(j + 

1 Ak—l,n 1 l,n 1 q Ak—l,n Ak—l,n—l 

+ /1 i-2,j + /1 i,j-2 + Z/1 i-l,j-l ~~ ^i-lj-l 

fe n k 

/ , ^10 + ^01 ,) -^tj l/MO + ^01 ) A ij 



1=1 m=l l=m 



(75) 



The A\™ yield the desired Rayleigh-Schrodinger coefficients E kn via the simple 
formula 

^n = -(-l) fc (< + <) • (76) 

These can be determined recursively via ([75]). The recursion must be initial- 
ized with 

A q = 5ko 5 n o , (77) 
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and solved at increasing k = 0, 1, 2, . . .; n = 0,1,2, ... ,k and, for each set 
k and n, with decreasing % = 2k — n, . . . , 0; j = 2k — n, ... , (omitting 
i=j=0). The procedure is most easily performed with the help of an algebraic 
computer program. The list of the first Rayleigh-Schrodinger coefficients up 
to k = 12 (n = 0, . . . , k) is given Table |l|. 



4.2 Large-order coefficients and resummation 



Working with Langer's formulation |2(| (which is related to Lipatov's |2T[ by 
a dispersion relation) and making use of known results for the isotropic an- 
harmonic oscillator, we shall derive the large-order behavior of perturbation 
expansion for the ground-state energy. 

The method is based on the path-integral representation of the quantum 
partition function 

Z = J VxVyexp[-A(x,y)} ^ exp(-(3E) (78) 

where 

A = / + ^ 2 dr [\ (i 2 + y 2 )+\ (x 2 + y 2 ) + I [x 4 + 2(1 - 6)xV + V 

~" 2 (79) 
is the Euclidean action corresponding to the Hamiltonian fl66|) . For g > 0, the 
system is stable and Z is real. On the other hand, if the coupling constant g is 
negative the system becomes unstable and Z develops an exponentially small 
imaginary part related to the decay-rate T of the ground-state resonance. 
The imaginary part of the ground-state energy may be obtained by taking 
the large j3 limit in ([78]), 

T „ 1 HmZ 

In the above equation the fact was used that ImZ oc exp(— (3) exp [— l/(a|g|)] 
is much smaller than ReZ = exp{— (3[1 + O(g)]}. For small g < 0, this 
imaginary part can be computed perturbatively in the anisotropic parameter 
S by the expansion around the isotropic instanton solution r c (r): 

( x \ — ( cosl P 



, i ,(r c + + { )V,r c = J — —- -.(81) 

y J \ smyj J 1^ cos^ y y \g\ cosh(r - r ) 
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k 


n 


Ekn 








1 


-1 

1 


rt 



o 

z 


1 




-1 /4 


2 





-9 


rt 

2 


1 


9/4 


2 


2 


-3/16 


3 





89 


3 


1 


-267/8 


3 


2 


177/32 


3 


3 


-11/32 


4 





-5013/4 


4 


1 


5013/8 


4 


2 


-9943/64 


4 


3 


2465/128 


4 


4 


-973/1024 


5 





88251/4 


5 


1 


-441255/32 


5 


2 


874757/192 


5 


3 


-216751/256 


5 


4 


171049/2048 


5 


5 


-20987/6144 


6 





-3662169/8 


6 


1 


10986507/32 


6 


2 


-327063703/2304 


6 


3 


81133049/2304 


6 


4 


-64093757/12288 


6 


5 


31487347/73728 


6 


6 


-4401593/294912 


7 





86716929/8 


7 


1 


-607018503/64 


7 


2 


32603176343/6912 


7 


3 


-40534191905/27648 


7 


4 


32089547489/110592 


7 


5 


-15794879119/442368 


7 


6 


4423646695/1769472 


7 


7 


-135064261/1769472 



8 





-18380724429/64 


8 


1 


18380724429/64 


8 


2 


-6932983533833/41472 


8 


3 


216189163547/3456 


8 


4 


-6865860756773/442368 


8 


5 


847050762955/331776 


8 


6 


-951145969207/3538944 


8 


7 


116411434099/7077888 


8 


8 


-151575359341 /339738624 


9 





537798950495/64 


9 


1 


-4840190554455/512 


9 


2 


785448510795415/124416 


9 


3 


-172109470699495/62208 


9 


4 


5484677894663731/6635520 


9 


5 


-2714832036203789/15925248 


9 


6 


1910496739715441/79626240 


9 


7 


-468820318449871 /212336640 


9 


8 


1223678377567247/10192158720 


9 


9 


-29878788733243/10192158720 


10 





-34427971992123/128 


10 


1 


172139859960615/512 


10 


2 


-757445337006448801/2985984 


10 


3 


95243865818145949/746496 


10 


4 


-26657139955813121627/597196800 


10 


5 


6619289843855618939/597196800 


10 


6 


-18688108386867852767/9555148800 


10 


7 


287396519063579707/1194393600 


10 


8 


-6016110774357344761/305764761600 


10 


9 


588950135128273907/611529523200 


10 


10 


-52319976745196951/2446118092800 
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11 





1196938085820951/128 


11 


1 


-13166318944030461/1024 


11 


2 


484953641311740249799/44789760 


11 


3 


-122498739278392549037/19906560 


11 


4 


273164737489274832749 / 1 10592000 


11 


5 


-8576021167229490768493/11943936000 


1 1 

1 1 


fi 

u 


4oo i 4oo i ooouzoy i uyooo/ oioou^youu 


11 


7 


-11098068163230668731/471859200 


11 


8 


786108601809348099491/305764761600 


11 


9 


-385711575108432009551/2038431744000 


11 


10 


7631665291905150913/905969664000 


11 


11 


-12593952190067271863/73383542784000 


12 





-179761724871375777/512 


12 


1 


539285174614127331 /1024 


12 


2 


-658487704407131831592119/1343692800 


12 


3 


83537029566207575386361/268738560 


12 


4 


-1870114571495468628478319/13271040000 


12 


5 


4208267075207850881725247/89579520000 


12 


6 


-16737064308714173333404777/1433272320000 


12 


7 


38348532616813953055927/17694720000 


12 


8 


-27234664511494120875149389/91729428480000 


12 


9 


1339555446357501564974269/45864714240000 


12 


10 


-53129831724844951147579/27179089920000 


12 


11 


70291727826145874647867/880602513408000 


12 


12 


-52920213881686076606297/35224100536320000 



Table 1: Coefficients E^ n in the perturbation series flT0| ) for the ground-state 
energy up to k = 12 (n = 0,. . . ,k). 
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For simplicity, we shall set r = in the sequel. In ( [81] ) we have separated 
out the rotation angle ip of the isotropic instanton in the (x, y)-plane and 
£, t] are the degrees of freedom, orthogonal to this rotation. Inserting the 
expansion (|81|) into the action fl79|) we obtain the expression 



.4 



4 5 2sin 2 (2<^) 




6 



cosh T 



dr A 



+ 1 



cosh T 



(82) 



where we have splitted of the action into the terms responsible for the leading 



contributions in an expansion of the form (|53| ) and a remainder O 
Then the 5-dependence of the transversal quadratic fluctuations belongs to 
the omitted terms. Expanding flT8|) in 5 and integrating out the quadratic 
fluctuations we obtain 



2- 



n=0 



2sin 2 (2<^) 



4 



exp 



3M 



where the angle integral can be done with 



2tt 



dip 2sin 2 (2(^) 



r 2 (n + i) 
r(2n+l) 



8 n 25 (n + -,n 



[1 + 0{g)\ 
(83) 

(84) 



The contribution and from the quadratic longitudinal and transversal 
fluctuations coincide with those appearing in the isotropic oscillator problem 
and are therefore known. With the isotropic classical action Aq c = 4/(3|g|) 
the well-known results are 



i 



2tt 







det' \-d 2 /dr 2 + 1 - 6/coshV) 



det(-d 2 /dr 2 + 1) 



-1/2 



I A 
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and 



fv 



1 3 A, 



Oc 



2tt 



det'(-d 2 /dr 2 + 1 - 2/coshV) 
det(-d 2 /dr 2 + 1) 



-1/2 



2^exp(-/3/2) (85) 



where we have used the partition function of the harmonic oscillator 

Z osc ee det(-d 2 /dr 2 + l)^ 2 = , 1 ^ «p(-/9/2) (86) 

2 smh(p/2) 

to normalize the determinants. In the upper determinants, the zero eigen- 
values are excluded. This fact is recorded by the prime. 

The longitudinal fluctuations £ contain a negative eigenmode, this be- 
ing responsible for the factor —i/2 and the absolut value sign, and a zero 
eigenmode associated with the translation invariance which is spontaneously 
broken by the special choice r = 0. The separation of this zero eigenmode 



in the framework of collective coordinates yields the factor (3 J Aq c /{2 , k) (see 
Chapter 17 in Ref. [0). Collecting the contributions of the negative and all 
positive eigenmodes one obtains the remaining factor \[V2. 

In contrast to the longitudinal case the transversal fluctuations r\ do not 
contain any negative mode. The transversal fluctuation operator has one 
zero eigenvalue due to the rotational invariance in the limit 5^0. The as- 
sociated eigenmode is extracted from the integration measure via the change 
of variables (Ell). The Jacobian of this coordinate treansformation can be de- 



duced from the isotropic system. It contributes the factor y3A 0c /(2n). The 
remaining factor 2 results from all other modes with positive eigenvalues. 

Collecting all contributions to the imaginary part of the ground-state 
energy (|80|), an cancellation of all (3 dependent factors leads to: 



ImE 9 -A - Y r—B [n+-,n + - 

/3exp(-/3)if n\ \ ' 

6 ~ (-25) n n / 1 1 

= - V t—B in + -,n + - 

irto n\ V 2' 2 

Finally, by means of the dispersion relation (153) we find the corresponding 
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large-order behavior of the coefficients in the expansion ([701): 

After having derived the large-order behavior of E^ n and the low-order per- 
turbation coefficients via the Bender and Wu-like recursions ([75]) and (|76|), 
we are in the position to resumme the g-series accompanying each power 5 n 
in the expansion (|7Cf). 

The remaining strong-coupling expansion follows from Symanzik scaling 

oo 

E(g,5)=J2^)9 {1 - 2m)/3 (89) 

m=0 

i. e., the power behavior in the strong- coupling limit is given by 

E(g,6) 9 ^ K (5)g l '\ (90) 

where the 1/3 coincides with that appearing in the one-dimensional oscillator 
problem. The (^-dependence enters by the prefactor k q . 

Combining (|17]) and the formulas (|60"D , ( |61~]) and (^), the resummation 
procedure must be worked through with the parameters 

i / \ 3 
b {n) = n + 



2 



3 



a 



4 

a = i. (91) 

In Figures [7| and || we have plotted the 5-dependence of the resummed ground- 
state energy E for two different values of the coupling constant g/4 and for 
various orders N. For a reference plot we have used the very accurate dotted 
curve which we have obtained numerically from the variational perturbation 
theory described in the next subsection. 



4.3 Variational perturbation theory 

It is useful to compare the above results with those of another recently- 
developed resummation procedure known as variational perturbation theory 
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Figure 7: Ground-state energy E of the anisotropic anharmonic oscillator 
with g/A = 0.1 as a function of the anisotropy 5. Shown are the resummed 
perturbation series for various order of approximation iV and the approxima- 
tion ^(Os) (see the text after Eq. (|99"D ) from the variational perturbation 
theory (VPT). Differences between W 5 (Q 5 ) and the exact ground-state en- 
ergy are existent only on a finer energy scale. 
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[for an introduction see Ref. [IT], Chapter 5]. Consider first the Rayleigh- 
Schrodinger expansion of the ground-state energy: 



E{g,5)=uf^Y,E lrn {25T( 9 -tf\ , (92) 

1=0 m=0 V U J 

where the Rayleigh-Schrodinger coefficients E\ m are obtained from the recur- 
sion relation (ff5|) via (|76|). 

The variation is done as follows: First, the potential is separated into an 
arbitrary harmonic term and a remainder: 

In contrast to ordinary perturbation theory, an interacting potential V- int is 
defined by 

V(x,y) = ^-(x 2 + y 2 )+V int (x,y) (94) 

and setting 

V int (x, y) = 9 - (pr 2 + r 4 - 28x 2 y 2 ) ; p = - (u 2 - fi 2 ) . (95) 
A perturbation expansion is now found in powers of g at fixed p and 8: 



k 



E k (g,8,p) = n£e l (p,S)\yp) . (96) 

The calculation of the new coefficients ei(p,8) up to a specific order k does 
not require much additional work, since they are easily obtained from the 
ordinary perturbation series (|9~2T) . We simply replace u by the identical ex- 
pression 

u = Vfi 2 + u 2 - Q 2 = y^ 2 + gp/2 , (97) 

reexpand E(g,8) in powers of g, and truncate the series after an order I > k. 
This yields the reexpansion coefficients 

ei(p, 8) = £ t E UW ( (1 7^ /2 ) (2 P n) l -i . (98) 

j=0n=0 \ J J 
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The truncated power series 



W k {g, 5, 0) := E k [g, 5, 2 (uj 2 - ft 2 ) / g] (99) 

is certainly independent of Q for k going to infinity. However, at any finite 
order it depends on Q. The optimal value of Q is found by calculating all 
extrema and the turning points. The smallest among these order-dependent 
points is used as an optimal trying value and is denoted by ftk(g, 6) . The asso- 
ciated energy Wk[g,5,£lk(g,5)] constitutes the desired approximation to the 
ground-state energy. In Figures [| and we have plotted the ^-dependence of 
Ws e f° r various anisotropy parameters 5 at the coupling constant g/A = 0.1. 
The shape of the curves depends little on 5. Only in the case of odd k does 
a minimum exist. For even k, there is no extremum and the optimal f2-value 
lies at a turning point. 

For an isotropic gx 4 -model, the precision of the variational perturbation 
method has been illustrated by a comparison with accurate numerical ener- 
gies |TJj . With increasing k the approach of H4 to the exact energy is quite 
rapid and its mechanism is well understood. 

In Table |2] we display the ground-state energies for odd k, which we have 
obtained for the anisotropic model at various 5 and g/A. The convergence to 
fixed energy values is comparable to the case of the simple g:r 4 -interaction. So 
we assume that these numbers coincide with the exact ground-state energy 
values at least up to the first four digits. 



5 Summary 

With the help of a simple model integral containing a quadratic and two 
quartic terms of different symmetry, we have investigated in detail the large- 
order behavior of the ^-dependent ^-series in a function f(g, 5) = J2k fk{8)g k 
for the region near the isotropic limit 5 — > 0. We have shown that the 
large-order behavior of fk{S) undergoes a crossover from the anisotropic to 
an isotropic regime near the order of perturbation theory k cross ~ l/\5\. 

In quantum mechanics, the extreme large-order behavior of perturbation 
theory for the anisotropic regime k\S\ 3> 1 is identical with earlier results of 
BBW || and Janke 0. In displaying the crossover-behavior we have gone 
beyond these earlier works. 
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Figure 9: The function W§ at constant coupling strength g/4 = 0.1 for 
various anisotropy parameters 5. 
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Figure 10: The function W 6 for g/A = 0.1 and various 5. 
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0/4 = 0.1 


k \S 


-2.5 


-1.5 


-0.5 


0.5 


1.5 


1 


1.222923 


1.19626 


1.167751 


1.137 


1.103438 


3 


1.217193 


1.192062 


1.164807 


1.134739 


1.100658 


5 


1.217109 


1.192032 


1.164801 


1.134734 


1.100607 


7 


1.217107 


1.192033 


1.164803 


1.134735 


1.100604 


9 


1.217107 


1.192034 


1.16481 


1.134736 


1.100604 


11 


1.217107 


1.192035 


1.16481 


1.134739 


1.100604 


g/A = 1.0 


k \5 


-2.5 


-1.5 


-0.5 


0.5 


1.5 


1 


1.969986 


1.88556 


1.791636 


1.684863 


1.559412 


3 


1.941934 


1.863112 


1.773978 


1.669261 


1.536823 


5 


1.941196 


1.862803 


1.773867 


1.669156 


1.535609 


7 


1.941172 


1.862806 


1.773888 


1.669172 


1.535454 


9 


1.941172 


1.862815 


1.773909 


1.669188 


1.535425 


11 


1.94118 


1.862823 


1.773924 


1.669199 


1.535418 



Table 2: Convergence of the ground-state energy in the variational pertur- 
bation expansion for various anisotropy parameters 5. 
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In particular, our resummation algorithm is shown to work very well 
in the vicinity of 8 = and for 5 > 0, the latter being relevant to the 
question of a stable cubic fixed point in field theory. With increasing coupling 
constant g/4, the error of the result for the ground-state energy becomes 
larger. However, for N = 6 (this is the largest available order for the (3- 
functions in quantum field theory, see Ref. ) and in the wide region 
5 G (—0.5,2) and g/4 G (0,1) the error remains smaller than 0.8%. The 
increasing error for large negative values of 5 can intuitively be understood 
by comparing the first two terms in the action fl32|): For 5 < 0, the "tunneling- 
paths" of extremal action are obviously straight lines along the two diagonales 
in the (x, y)-plane (if = tt/4). Along these diagonal rays, the basic factor 
exp[— l/(c|g|)] related to the decay-rate disappears for 5 — > —2, and the 
ensuing expansion of (|78|) in powers S n becomes meaningless. An improved 
fit for 5 < can be obtained by choosing larger values of the large-order 
parameter a. In Figures [11] and [12| we display the result for o = 3 and 
N = 6, where for g/4 = 0.1 the accurate and the resummed curve coincide. 

To obtain the correct description of the neighbourhood of the isotropic 
system 5 = we have used the method developed in the context of an 
anisotropic quantum field theory in 0: By replacing the series J2k fk(8)g k by 

Sfc fkn 9 k 8 n and resumming the g-series accompanying each power S n , we 
obtain very good results for the model integral and the ground-state energy 
of the anisotropic anharmonic oscillator. In this way our results justify the 
earlier field theoretic analysis and should be useful for understanding similar 
problems in other systems. 
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E 




8 



Figure 11: The same functions as in Fig. [7], but with the large-order 
parameter a = 3 (explained in the text). 
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Figure 12: The same functions as in Fig. |8|, but with the large-order 
parameter a = 3. 
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